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Abstract 

Objectives: It is critical to implement effective multiple countermeasures to 
mitigate or retain the spread of pandemic influenza. We propose a mathematical 
pandemic influenza model to assess the effectiveness of multiple countermea- 
sures implemented in 2009. 

Methods: Age-specific parameters, including the transmission rate, the propor- 
tion of asymptomatic individuals, the vaccination rate, the social distancing rate, 
and the antiviral treatment rate are estimated using the least-square method 
calibrated to the incidence data. 

Results: The multiple interventions (intensive vaccination, social distancing, 
antivrial treatment) were successfully implemented resulting in the dramatic 
reduction in the total number of incidence. 

Conclusion: The model output is sensitive to age-specific parameters and this 
leads to the fact that a more elaborate age group model should be developed and 
extensive further studies must be followed. 



1. Introduction 

Influenza imposes serious social and economic 
burden to many countries all around the world [1]. In the 
US, seasonal influenza results in 200,000 hospitaliza- 
tions and 36,000 deaths annually, thus its economic 
burden amounts up to $87.1 billion per year [2]. In 
France, the economic loss due to the absence from 
workplace caused by influenza is approximately 13,076 



French Francs (about €2,431) in a year [3]. People in 
UK reportedly miss on average 2.8 workdays because of 
influenza [4]. In Germany, the per unit cost of an 
influenza case in 1996—97 was 1,777 Deutsche Mark 
(about €1,105.63) [5]. According to Szuch [6], the 
productivity loss per unit due to missing workdays 
because of influenza was to range from €1,379 to 
€6,991 and from €482 to €1,409 due to direct infec- 
tion. In addition to direct economic burden, spread of 
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influenza can cause psychological burden that is not 
reported in the statistical records on economic loss. 

Therefore, government and public health officials in 
many countries have made their efforts to resist against 
the spread of influenza, which is especially the case 
when it comes to pandemic influenza. The counter- 
measure strategy includes vaccination, social distancing, 
and anti-viral treatment. Given the limited amount of 
available resources, it is critical to find the most effec- 
tive strategy or multiple strategies before the influenza 
takes place. Moreover, it is crucial to assess the effec- 
tiveness of these countermeasures afterwards since it 
would provide invaluable information for the future 
influenza plan. Mathematical modeling is useful for both 
aims. Using mathematical models, we can simulate how 
the epidemic would change when we utilize specific 
countermeasures. Also, we can calibrate it using 
empirical data and assess the effectiveness of counter- 
measures which was implemented in the past. The latter 
approach is especially conducive when we have only 
data which includes the impact of a variety of counter- 
measures [7,8]. Assessing the effectiveness of each 
countermeasure would increase the possibility that we 
can handle the influenza more efficiently for the future. 
It was the case of SARS in 2003 where models were 
built based on past data and appropriate intervention 
strategies were implemented based on the predictions 
that the models produced. 

This stud focuses on the case of 2009 Influenza A 
(H1N1) in the Republic of Korea (hereafter Korea). 
Influenza A (H1N1), a mutant of swine flu which is 
known to appear first in Mexico in 2009 and spread to 
the whole world, has been a serious public health 
problem as well as social and economic ones throughout 
the globe [9,10]. In the US, according to CDC, about 
600 million people, which amounted up to 20% of total 
US population, were reportedly to be infected by the 
influenza. In Korea, after a traveler to Mexico was 
identified to be infected in April 2009, the number of 
infected was peaked in November. The Korean health 
authorities implemented a vaccination program to the 
hospital personnel from October 27 and expanded the 
coverage of vaccination to the general public from 
November 11 [11]. This intervention turned the diffu- 
sion trend downward, and the Influenza A (H1N1) was 
finally declared to be eliminated from Korea in October 
2010. The peculiarity of this disease was the high in- 
fectious rate of the younger age group and low rate of 
the older (65 and over) age group [12]. It is believed that 
the older age group get partially immune when they 
have experienced Spanish Influenza in the past [13]. 

As the case above and others show, when it comes to 
assessing the effectiveness of countermeasures, the age 
structure of population should be taken into consider- 
ation. It is no wonder because people in different age 
groups can be justifiably assumed to have different 
health conditions and different contact rates which come 
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Figure 1. Contact matrix between and within age groups. 



from different social and economic behaviors. There 
have been many previous researches about the effec- 
tiveness of countermeasures include age structures into 
their models [2,14,15,16,17]. 

This study presents a mathematical model with three 
age groups of the pandemic of Influenza A (H1N1) in 
2009. Also, using the incidence data in Korea, we carry 
out parameter estimations where the best-fitted param- 
eters are sought by the least- square method. The effec- 
tiveness of three intervention strategies, which are age- 
specific vaccination, social distancing, and antiviral 
treatment, is compared by calculating the basic repro- 
duction number, R 0 . 

2. Materials and methods 

2.1. Influenza pandemic transmission model 
with age groups 

We integrated the age structure of the Korean popu- 
lation to the influenza transmission model, based on the 
2009 Census data [18]. The Korean population was then 
divided into the following three age groups: Group 1, 
0—19 years; Group 2, 20—64 years; and Group 3, >65 
years. Further, each age group (indexed by i) is classified 
into eight epidemiological states, namely, susceptibility 
(Si), effectively vaccinated but not yet protected (V t ) 9 
latent (E t ), symptomatic and infectious (I f ), asymptom- 
atic and infectious (A t ), hospitalized (J t ), recovered (Rf), 
and dead (D t ). Susceptible individuals in age group i are 
exposed to the influenza virus at the force of infection: 

where ft is the transmission rate of age group i, which is 
assumed to be constant within age groups. The total 
population size N(t) is given by: 
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N(t)=J2 s k(t) + V k {t) +E k {t) +4(0 +A k (t) + J k (t) 
k=\ 

The force of infection consists of the contact rates 
0 z y, which are the age-specific contact rates modeled 
based on a study describing self-reported age-specific 
contact rates for the spread of respiratory infections 
[14]. The contact rate matrix is highly assortative with 
higher mixing rates within each age group than be- 
tween groups. Contact rates among 20—64 years old 
are the highest and rates among seniors (>65 years 
old) are the lowest (see Figure 1 and the contact matrix 
(2)). 



the epidemic onset with a vaccination rate v(f). That is, 
v(t) = 0, whenever t < t*. Age-specific vaccine efficacy 
is denoted by a t . Successfully vaccinated individuals 
progress to be protected while ineffectively vaccinated 
individuals remain susceptible to infection. Vaccinated 
but not yet protected individuals (V f ) may still be 
infected with influenza at the age-dependent force of 
infection A f as described above. 

Age-specific parameter values are described in 
Tables 1 and 2. The population is assumed to be 
completely susceptible at the beginning of the 
epidemic. The system of differential equations 
that describes our influenza transmission model is 
given by: 



St (0 = - £(0ft(0 E [4>ij{bAj{t) + (1 - u^Ijif)}} /N(t) - v f (0£(0 

7=1 

Vi = Vi{t)Si{t) - (1 - iTiMtMt) E [</>ij{bAj + (1 - Ui )Ij{t)}]/N{t) 

7=1 

< £ f = {£(0 + (1 " °i)Vt{t)} E h(t)/N(t) - kEi{t) (1) 

7=1 

I i {t)=kpE i {t) - + y + ti)Ii(t) 
A i (t)=k(\-p)E i (t)-yA i (i) 
j i {t) = a i I i {t)-{6 i + d i )J i {t) 
Ri(t) = yAi(t) + (y + ti)Ii(t) + SiJiit) 

Vb i (t)=5 i j i (t) 



Latent individuals E t progress to the infectious class 
Ii at the rate k {Ilk is the mean latent period). Infectious 
individuals are hospitalized at the age-specific mean 
rates a t and recover at the mean rate y. Hospitalized 
individuals either recover at the constant rate 6 t or die 
from influenza at the age-specific rate 5*. While the age- 
specific hospitalization rates are adjusted using esti- 
mates of the probability of hospitalization given clinical 
illness by age group, the recovery rate # z is assumed to 
be constant across all age groups for simplicity. 
Recovered individuals are assumed to remain protected 
for the duration of the epidemic. Vaccination, social 
distancing, and antiviral treatment are implemented after 
the incidence reaches the peak. For instance, vaccination 
is administered to susceptible individuals t* days after 



<i>ir 



131.8 39.2 6.4 
39.2 268.9 34.8 
6.4 34.8 76.0 



*J=1,2,3 



(2) 



Mii= m*»(^ + HLzPl) iJ=w (3) 



N \ai + y y 
_NM iJ (p{\-u i ) b(l-p) 



N \ai + v + U 



y 



U=W (4) 



The basic reproduction number, in the absence of 
interventions, % is given by the maximum eigenvalue 
of Equation (3) and similarly, the controlled basic 
reproduction number, 3? c can be computed from the 
Equation (4). 



Table 1. Age-specific parameter values used in numerical simulation 



Age group 


Peak incidence 


Total incidence 


Hospitalization rate 


Mortality rate 


Vaccine efficacy 


0-19 yr 


115,777 


552,802 


1.19 


0.007 


0.8 


20-59 yr 


25,778 


198,878 


0.77 


0.049 


0.8 


>60 yr 


1516 


12,081 


4.06 


0.483 


0.6 
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Table 2. Parameter values 



Parameters 


Description 




Value 






Refs 


T 


Recovery rate for 
infectious individuals 


7/4 








[1] 


e 


Recovery rate for 
hospitalized individuals 


2.38 








[1] 


k 


Rate of progression from 
latent to infectious 
individuals 


7/(1.2) 








[1] 


b 


Relative infectiousness of 
asymptomatic cases 
compared with infectious 
cases 


0.142 








[1] 


ft 


Probability of transmission 


0.0602 


2=1 






Data fitted 




per contact 
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0 




0.12 


2=1 
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0 1 

U. 1 


7* — 3 
6 — J 




Ui 


Social distancing rate 


0.03 


2=1 


0.15 


2=1 


Data fitted 






0.05 


2 = 2 


0.1 


2 = 2 








0 


2 = 3 


0.1 


2 = 3 




ti 


Antiviral treatment rate for 
hospitalized individuals 


0.1 


2=1 


0.25 


2=1 


Data fitted 






0.15 


2 = 2 


0.2 


2 = 2 








0.1 


2 = 3 


0.2 


2 = 3 





3. Results 

3.1. Parameter estimation 

Using the parameter values estimated through the 
least-squares method (Tables 1 and 2), the model output 
is illustrated in Figure 2. In our model, the infected 
incidence peak number is 1 15,780 (Group 1) at 44 week. 
In addition, the number for Group 2 is 25,767 at 45 week 
and for Group 3 it is 1514 at 44 week (Table 1). Relative 
error ( {actual val ™)-W™ lated values A wlt h data and simula- 

y actual value 1 

tion number is 0.00003 (Group 1), 0.0004 (Group 2), and 
0.0013(Group 3). 

Figure 3 is described simulation of non-age-groups 
with estimated parameters. In order to compared with 
total infected incidence number data per month with our 
simulation per week, we calculated (about every 4 
weeks) the cumulative infection number per week in our 
model (1). The value of bootstrap method 95% confi- 
dence interval was calculated and displayed by resam- 
pling 100,000 records. Bootstrap Method determines 
how accurate our estimation value by the number of 
times random resampling. Figure 3 shows that influenza 



data are almost in the 95% confidence interval, so it is 
explained that our model is reasonable by results of 
Figure 2 and Figure 3. 
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Figure 2. Age-specific incidence data of 2009 H1N1 influ- 
enza (bar graph) and its best-fitted simulation results (curves). 
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Figure 3. Total incidence data (*) and 95% confidence in- 
terval (gray area) using the bootstrap method. 



3.2. Effectiveness of intervention strategies 

Figure 2 presents the result of simulation to find the 
parameters of our model to fit the data. Then, using this 
calibrated model, we can evaluate the effectiveness of 
each control measure (vaccination, social distancing, 
and antiviral treatment). Parameter values used in 
simulation are shown in Table 2; parameter values of 
vaccination are 0.02 (Groups 1 and 3) and 0.03 (Group 2), 
the parameter value of social distancing is 0.02, and that 
of antiviral treatment is 0.2 since 46 week. 

The effectiveness of intervention strategies is illus- 
trated in Figures 4 and 5. Figure 4 compares the 
epidemic curves with and without control for each age 
group where solid curves are the results without controls 
and dotted curves are the results with controls. The peak 
of the solid curve (without control) is higher and earlier 
than the peak of the dotted curve (with control). 




18 23 27 31 36 40 45 49 1 6 10 14 18 23 27 32 36 
Time (wk) 



Figure 4. Age-specific incidence with controls (dotted curves) and without control (solid curves). 



x 10 4 Vaccine only x -|g 4 Social distancing only x -jq 4 Treatment only 




Time (Mo) 

Figure 5. Age-specific incidence is illustrated under three countermeasures; with controls (broken curves) and without controls 
(unbroken curves). The results under (A) vaccination, (B) social distancing, and (C) antiviral treatment. 
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Figure 5 demonstrates the simulation results when 
only one intervention strategy is implemented. 
Figure 5A shows the vaccination-only strategy, 
Figure 5B for social distancing only, and Figure 5C is 
for antiviral treatment only. It needs to be noted here 
that these interventions are implemented two months 
before the peak, because the effect is not appeared when 
controls are given one month before the peak. In the 
vaccination only case, the vaccination level for Group 1 
is 0.12 and the ones for other groups are 0.1. The cu- 
mulative infected incidence numbers in this case are 
1,180,400 (Group 1), 276,400 (Group 2), and 18,200 
(Group 3). In the social distancing only case 
(Figure 5B), the cumulative infected incidence numbers 
are 1,452,400 (Group 1), 498,000 (Group 2), and 36,800 
(Group 3). In addition, in the antiviral treatment only 
case, the numbers are 1,410,220 (Group 1), 475,700 
(Group 2), and 31,800 (Group 3). As can be seen in 
Figure 5, the infected incidence number is most reduced 
when the vaccination only strategy was implemented, 
the next effective was social distancing only, and the 
antiviral treatment only was the least effective. 

4. Discussion 

Devising effective countermeasures against influenza 
is one of the major concerns in public health officials. It 
is not only for economic cost-effectiveness but also for 
psychological stableness of people in society. Since it is 
impossible to conduct experiments on the spread of 
influenza in real-world settings, mathematical models 
are of great use to tackle the issue. Its usefulness 
stretches not only to forward-looking prediction of the 
future transmission under a variety of conditions but 
also to backward-looking assessment of combined anti- 
influenza measures in the past. 

This study presented a mathematical model of 
Influenza A (H1N1) with age three groups. Parameter 
estimation is carried out using the least square method to 
the 2009 pandemic influenza incidence data in Korea. 
The relative errors of incidence peak were 0.00003, 
0.0004, and 0.0013 for each group, and the fit of model 
with the 95% of confidence which was calculated by the 
Bootstrap method was given. The basic reproduction 
number R0 was 1.44 and it reduced to 1.14 after 
implementation of intensive interventions. 

It was also explored how much influence each control 
has on the infected incidence. The most successful 
intervention strategy was vaccination and social 
distancing followed by antivrial treatment when each of 
intervention strategy is implemented separately. The 
most effective intervention would be a mixed strategy 
which combines vaccination, social distancing and 
antiviral treatment all together. 

It is suggested that more elaborate age groups need to 
be incorporated in the model for our future study. 



Although the present study utilized three age groups, 
age groups can be more finely subdivided based on a 
health or social condition. For example, people in their 
20s and 50s can be justifiably assumed to have different 
resistance to influenza. Also, they may have different 
ways of social interactions and these may have impact 
on their possibility to be infected. Reflecting this factors 
and utilizing more elaborate age groups will improve the 
plausibility of the model and enable us to more accu- 
rately assess the effectiveness of countermeasures. 
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